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ABSTRACT 

This paper presents the fundamental properties 
and some numerical results of the solution of the 
diffusion equation of an impulsive cosmic-ray point 
source in an uniform, unbounded and spherically 
symmetrical moving medium. 


I. INTRODUCTION The diffusion-convection(D-C^ function is 
an elementary composite function of the solution of the D-C 
equation for the particles injected impulsively from a diffu- 
sive point source into a uniform, unbounded and spherically 
symmetrical medium. It is the analytic solution derived by 
the dimensional method for the propagation equation of solar t(J 
cosmic rays in the heliosphere, i.e. the interplanetary space. 
Because of the introduction of convection effect of solar 
wind, a nonhomogeneous term appears in the propagation equa- 
tion, it is difficult to express its solution in terms of 
the ordinary special functions. The research made so far has 
led to a solution containing only the first order approxi- 
mation of the convection effect. Undoubtedly it is insuffi- 
cient. to study only the propagation of the particle with not 
too high energy. The solution we get with the dimensional 
method, up to date, is the most general solution of propa- 
gation equation in a uniform, unbounded and spherically sym- 
metrical medium. It includes the higher orders of approxima- 
tion for the convection effect, and has been used in discus- 
sing various kinds of propagation effects of solar cosmic 
ray, making propagation corrections and evaluating the equi- 
valent diffusion coef ficients.** 3 The solution rnay also have 
its value for reference in the discussion of the diffusion 
in che ordinary moving medium. 

It is nece33ary to point out that this solution is 
valid only in the case of a constant K. For the low energy 
particles, the energy of particles and also their diffusion 
coefficients cnange as a result of the adiabatic expansion 
lo 3 s and other processes. In tnis case, such limitations 
can be retrieved by the equivalenb diffusion coefficient 
reflecting the average nature of propagation space. 

The diffusion equation in a uniform, unbounded and 

spherically symmetrical medium is 
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whore U is tne number density of diffusing particle, C is a 
constant relating to tho energy spectrum index. In the space 
far from the source, the solution satisfying condition 
* Presently at Department of Physics, University of New 
Hampshire, Durham, NH 0382ii, U.S.A. 
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U(r, t=o) = -^p£(r-r,' and Eq. (1.1' is 
where 4> is •diffusion-convection function* which the paper 


(i-i) 


is going to discuss in detail. It can be expressed as a func- 
tion of two dimensionless parameters as follows: Diffusion 
parameter Z— r*/(l|iCt), Convection parameter for time 

X— (2C - 1)V( t/K) . Disregarding 4 > , solution (1.2) is like 
a kind of the diffusion caused by a source moving outward 
with the convective velocity V. The equation of is 

'*T7. <*->■$* Cl5 > 

In case of no convection, V = 0 and<f> = l, 3olution(1.2' is a 
function of a single parameter Z. When convective effect is 
not too strong, 4> can be expanded into a cower series of 
power series of parameter X: 7) ) — IX^F(1,Z', 

F(1,Z' is called the 1-th order influence function of 
vection, and its equation is 

i >■*[-» *4rj F'a.?) F(KO 

Seeing that the solution must be finite at the origin, 
expression of F|l,Z' w can be taken as 
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II. MAIN METHODS A(n,m,k) is a multiple series relating 
to m,and k, the lowest recurrence values for m and k are 
all zero. For the derivation of the general recurrence 
formula of A(n,m,k), it is instructive to evaluate first 
from the lower-order coefficient A(n,m-l,k' the next one 
<\(n,m,k' byrecurring k. Then, the general formula of A(n,m,k) 
can be yielded as: r -j 
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where coefficient A(n) is determined by the initial con- 
dition. 

The multiple series s(n,m,k' can be transformed into 
a series composed of the higher-order T* function (function 
and Cr) . 

From the initial condition at t * o and the asymptotic 
value (*T ) of function as Z-»« , the recurrence formula of 
coefficient A(l' can be yielded : 

«i»s’ 4c*)(4) H S<o,i)-o^ C*-») 
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in n h i G1)or ? I rder coefficient A is difficult to express 

late HoweVer > ifc is convenient to calcu- 

laze Ail with the recurrence formulae and ( z*5 ) bv 

uflnr« ter *^ I h Can v b8 shown that A(1) is an alternative seq- 
uence, and its absolute value decreases when 1 increases. 

i ' orraul a ( 1,5 ) is a nonhomogeneous equation. Its non- 
homogeneous terms are composed of the same functions of lower 

S’ formula of function F<1,Z) can be derived from 

formulae ( »-6 ) and ( 2-1 )j 

F(l.I) =£ o C***, * ) <.*- 4 ) 

The asymptotic value of f(n,m,Z) when t * o or Z-*~ i s 

-rc4> 2*^ ^f(g) (^) H so,,«) OT; 

Therefrom the recurrence formula ( 2-5 ^ of coefficient A(l> 
can be formulated. 


Il l , DIPFUSION-^OIWaCTIQN F U NCTION 

. .. i ' or convenience of discussion, function <f> can be re- 
written as <*> 

4 > = I Acn)ac«>*«) x"Y* -f ) (a -0 

where Y is another dimensionless parameter, namely, the con- 
vection parameter for space j „ . . Vr *h 

Y~ ( 2G - D-pft- = XZ^ 

Pig. 1. illustrates the value of n for calculation of 
* as a function of Y and T, when the relative magnitude of 
the last term os less than 10 . It can be seen that the 

term number of convergence increases with the value of Y 
and T. Functions is probably divergent as X > 3roo. So, we 
often restrict our discussion to propagation problem within 
the range X «s 3.oo. 


F.q.1. 


Abscissa is the dimensionless time T> ordinate 
is the number of terms summing up the function 
<p accurate to 10“\ 



iy. VARIATION OF DENSITY U Variation of n-c function 
with the convection parameter for space Y and the dimension- 
less time t( = Ij-Kt/r 1 ) have already been shown in Pig. 1. 
of Cal . Plotted here is the variation of diffusing particles 
density U with T under the different conditions of solar 
wind convecbions (Pig. 2.). A set of curves on the left 
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side are the results calculated by Formula ( 1*2. ), talcing 
2C - 1 = Ip, or the differential -momentum- spectrum index of 
particles to be 5.5. q uch case corresponds to that of the 
isotropic propagation. But in fact, propagation of cosmic 
rays particles in interplanetary space is anisotropic, the 
diffusion coefficient along the magnetic line YL f , is diffe- 
rent from the transverse one Kj. . In this ca3e, propaga- 


tion equation becomes 


u “ 


i* ^ (—)' «*? l- 2 t ,+ % i- W) 


where is the magnetic azimuth of observation point from 
the source, and H i3 the modulation parameter, Vr 

2KT 

Substitute k/, for K when Z is calculated. 

Solution (I4..I) of the anisotropic propagation equa- 
tion is drawn by a set of curves of the right side of Fig2 
when 6 =0*. It can be shown that the intensity of diffu- 
sing particles decreases and the maximum time moves up as 
the convection effect intensifies. Besides, the time scale 
of tne isotropic D-0 propagation is longer than that of 
the anisotropic one* -Hence, it is possible to form an over 
estimation of equivalent diffusion coefficient if obser- 
vational data is fitted merely by the former. However, bet- 
ter results can be got by using the model of anisotropic 



Dimcnsionle* lime T 

(a) Isotropic D-C (b) anisotropic P-C 
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